Photoemission Spectra from the Extended Koopman's Theorem, Revisited.

The Extended Koopman's Theorem (EKT) provides a straightforward way to compute charged excitations from any level of theory. In this work we make the link with the many-body effective energy theory (MEET) that we derived to calculate the spectral function, which is directly related to photoemission spectra. In particular, we show that at its lowest level of approximation the MEET removal and addition energies correspond to the so-called diagonal approximation of the EKT. Thanks to this link, the EKT and the MEET can benefit from mutual insight. In particular, one can readily extend the EKT to calculate the full spectral function, and choose a more optimal basis set for the MEET by solving the EKT secular equation. We illustrate these findings with the examples of the Hubbard dimer and bulk silicon.

A diagonalization-free optimization algorithm for solving Kohn-Sham equations of closed-shell molecules.

A local optimization algorithm for solving the Kohn-Sham equations is presented. It is based on a direct minimization of the energy functional under the equality constraints representing the Grassmann Manifold. The algorithm does not require an eigendecomposition, which may be advantageous in large-scale computations. It is optimized to reduce the number of Kohn-Sham matrix evaluations to one per iteration to be competitive with standard self-consistent field (SCF) approach accelerated by direct inversion of the iterative subspace (DIIS). Numerical experiments include a comparison of the algorithm with DIIS. A high reliability of the algorithm is observed in configurations where SCF iterations fail to converge or find a wrong solution corresponding to a stationary point different from the global minimum. The local optimization algorithm itself does not guarantee that the found minimum is global. However, a randomization of the initial approximation shows a convergence to the right minimum in the vast majority of cases.

Optical spectra of 2D monolayers from time-dependent density functional theory.

The optical spectra of two-dimensional (2D) periodic systems provide a challenge for time-dependent density-functional theory (TDDFT) because of the large excitonic effects in these materials. In this work we explore how accurately these spectra can be described within a pure Kohn-Sham time-dependent density-functional framework, i.e., a framework in which no theory beyond Kohn-Sham density-functional theory, such as GW, is required to correct the Kohn-Sham gap. To achieve this goal we adapted a recent approach we developed for the optical spectra of 3D systems [S. Cavo, J. A. Berger and P. Romaniello, Phys. Rev. B, 2020, 101, 115109] to those of 2D systems. Our approach relies on the link between the exchange-correlation kernel of TDDFT and the derivative discontinuity of ground-state density-functional theory, which guarantees a correct quasi-particle gap, and on a generalization of the polarization functional [J. A. Berger, Phys. Rev. Lett., 2015, 115, 137402], which describes the excitonic effects. We applied our approach to two prototypical 2D monolayers, h-BN and MoS2. We find that our protocol gives a qualitatively good description of the optical spectrum of h-BN, whereas improvements are needed for MoS2 to describe the intensity of the excitonic peaks.

Clifford Boundary Conditions: A Simple Direct-Sum Evaluation of Madelung Constants.

We propose a simple direct-sum method for the efficient evaluation of lattice sums in periodic solids. It consists of two main principles: (i) the creation of a supercell that has the topology of a Clifford torus, which is a flat, finite, and borderless manifold; (ii) the renormalization of the distance between two points on the Clifford torus by defining it as the Euclidean distance in the embedding space of the Clifford torus. Our approach does not require any integral transformations nor any renormalization of the charges. We illustrate our approach by applying it to the calculation of the Madelung constants of ionic crystals. We show that the convergence toward the system of infinite size is monotonic, which allows for a straightforward extrapolation of the Madelung constant. We are able to recover the Madelung constants with a remarkable accuracy, and at an almost negligible computational cost, i.e., a few seconds on a laptop computer.

Many-Body Effective Energy Theory: Photoemission at Strong Correlation.

In this work, we explore the performance of a recently derived many-body effective energy theory for the calculation of photoemission spectra in the regime of strong electron correlation. We apply the theory to paramagnetic MnO, FeO, CoO, and NiO, which are typical examples of strongly correlated materials and, therefore, a challenge for standard theories. We show that our method opens a correlation gap in all of the oxides studied without breaking the symmetry. Although the materials seem similar, we show that an analysis of the occupation numbers reveals that the nature of the gap is not the same for these materials. Overall, the results are very promising, although improvements are clearly required, since the band gap is overestimated for all of the systems studied. We indicate some possible strategies to further develop the theory.

Unphysical Discontinuities in GW Methods.

We report unphysical irregularities and discontinuities in some key experimentally measurable quantities computed within the GW approximation of many-body perturbation theory applied to molecular systems. In particular, we show that the solution obtained with partially self-consistent GW schemes depends on the algorithm one uses to self-consistently solve the quasiparticle (QP) equation. The main observation of the present study is that each branch of the self-energy is associated with a distinct QP solution and that each switch between solutions implies a significant discontinuity in the quasiparticle energy as a function of the internuclear distance. Moreover, we clearly observe "ripple" effects, i.e., when a discontinuity in one of the QP energies induces (smaller) discontinuities in the other QP energies. Going from one branch to another implies a transfer of weight between two solutions of the QP equation. The cases of occupied, virtual, and frontier orbitals are separately discussed on distinct diatomics. In particular, we show that multisolution behavior in frontier orbitals is more likely if the HOMO-LUMO gap is small.

Distributed Gaussian orbitals for the description of electrons in an external potential.

In this work, we demonstrate the viability of using distributed Gaussian orbitals as a basis set for the calculation of the properties of electrons subjected to an external potential. We validate our method by studying one-electron systems for which we can compare to exact analytical results. We highlight numerical aspects that require particular care when using a distributedGaussian basis set. In particular, we discuss the optimal choice for the distance between two neighboring Gaussian orbitals. Finally, we show how our approach can be applied to many-electron problems.

Green Functions and Self-Consistency: Insights From the Spherium Model.

We report an exhaustive study of the performance of different variants of Green function methods for the spherium model in which two electrons are confined to the surface of a sphere and interact via a genuine long-range Coulomb operator. We show that the spherium model provides a unique paradigm to study electronic correlation effects from the weakly correlated regime to the strongly correlated regime, since the mathematics are simple while the physics is rich. We compare perturbative GW, partially self-consistent GW and second-order Green function (GF2) methods for the computation of ionization potentials, electron affinities, energy gaps, correlation energies as well as singlet and triplet neutral excitations by solving the Bethe-Salpeter equation (BSE). We discuss the problem of self-screening in GW and show that it can be partially solved with a second-order screened exchange correction (SOSEX). We find that, in general, self-consistency deteriorates the results with respect to those obtained within perturbative approaches with a Hartree-Fock starting point. Finally, we unveil an important problem of partial self-consistency in GW: in the weakly correlated regime, it can produce artificial discontinuities in the self-energy caused by satellite resonances with large weights.

Signatures of Wigner localization in one-dimensional systems.

We propose a simple and efficient approach to study Wigner localization in one-dimensional systems using ab initio theory. In particular, we propose a suitable basis for the study of localization which consists of equally spaced overlapping gaussians. We illustrate our approach with full-configuration interaction which yields exact results for a given basis set. With our approach, we were able to study up to 8 electrons with full-configuration interaction. Finally, we propose the total-position spread tensor and the total electron entropy as convenient quantities to obtain signatures of Wigner localization.

Many-body perturbation theory and non-perturbative approaches: screened interaction as the key ingredient.

Many-body perturbation theory is often formulated in terms of an expansion in the dressed instead of the bare Green's function, and in the screened instead of the bare Coulomb interaction. However, screening can be calculated on different levels of approximation, and it is important to define what is the most appropriate choice. We explore this question by studying a zero-dimensional model (so called 'one-point model') that retains the structure of the full equations. We study both linear and non-linear response approximations to the screening. We find that an expansion in terms of the screening in the random phase approximation is the most promising way for an application in real systems. Moreover, by making use of the nonperturbative features of the Kadanoff-Baym equation for the one-body Green's function, we obtain an approximate solution in our model that is very promising, although its applicability to real systems has still to be explored.

Gauge-Invariant Formulation of Circular Dichroism.

Standard formulations of magnetic response properties, such as circular dichroism spectra, are plagued by gauge dependencies, which can lead to unphysical results. In this work, we present a general gauge-invariant and numerically efficient approach for the calculation of circular dichroism spectra from the current density. First we show that in this formulation the optical rotation tensor, the response function from which circular dichroism spectra can be obtained, is independent of the origin of the coordinate system. We then demonstrate that its trace is independent of the gauge origin of the vector potential. We also show how gauge invariance can be retained in practical calculations with finite basis sets. As an example, we explain how our method can be applied to time-dependent current-density-functional theory. Finally, we report gauge-invariant circular dichroism spectra obtained using the adiabatic local-density approximation. The circular dichroism spectra we thus obtain are in good agreement with experiment.

Fully Parameter-Free Calculation of Optical Spectra for Insulators, Semiconductors, and Metals from a Simple Polarization Functional.

We present a fully parameter-free density-functional approach for the accurate description of optical absorption spectra of insulators, semiconductors, and metals. We show that this can be achieved within time-dependent current-density-functional theory using a simple dynamical polarization functional. We derive this functional from physical principles that govern optical spectra. Our method is truly predictive because not a single parameter is used. In particular, we do not use an ad hoc material-dependent broadening parameter to compare theory to experiment as is usually done. Our approach is numerically efficient; the cost equals that of a calculation within the random-phase approximation.

Reduced density-matrix functional theory: Correlation and spectroscopy.

In this work, we explore the performance of approximations to electron correlation in reduced density-matrix functional theory (RDMFT) and of approximations to the observables calculated within this theory. Our analysis focuses on the calculation of total energies, occupation numbers, removal/addition energies, and spectral functions. We use the exactly solvable Hubbard dimer at 1/4 and 1/2 fillings as test systems. This allows us to analyze the underlying physics and to elucidate the origin of the observed trends. For comparison, we also report the results of the GW approximation, where the self-energy functional is approximated, but no further hypothesis is made concerning the approximations of the observables. In particular, we focus on the atomic limit, where the two sites of the dimer are pulled apart and electrons localize on either site with equal probability, unless a small perturbation is present: this is the regime of strong electron correlation. In this limit, using the Hubbard dimer at 1/2 filling with or without a spin-symmetry-broken ground state allows us to explore how degeneracies and spin-symmetry breaking are treated in RDMFT. We find that, within the used approximations, neither in RDMFT nor in GW, the signature of strong correlation is present, when looking at the removal/addition energies and spectral function from the spin-singlet ground state, whereas both give the exact result for the spin-symmetry broken case. Moreover, we show how the spectroscopic properties change from one spin structure to the other.

Gauge-invariant calculation of static and dynamical magnetic properties from the current density.

In this work we solve two problems related to the calculation of static and dynamical magnetic properties with ab initio theories. First, we show that the dependence of the dynamical magnetic dipole moment on the reference point of the multipole expansion and on the gauge origin of the vector potential have a clear physical significance. They are due to a dynamical electric dipole moment and an electric field, respectively. Both are fully determined by the experimental setup and do not pose any fundamental problem, contrary to what is commonly assumed. Second, in the static case, any dependence on the gauge origin is an artifact of the computational method. We show that the artificial dependence on the gauge origin can be removed in an elegant way by the introduction of a sum rule that puts the diamagnetic and paramagnetic contributions on equal footing. Our approach can be applied to calculate any magnetic observable that can be derived from the current density, and can be used in combination with any ab initio theory from which it can be obtained. To illustrate our method we apply it here to time-dependent current-density-functional theory for the calculation of static and dynamical magnetizabilities of molecules.

A rational reduction of CI expansions: combining localized molecular orbitals and selected charge excitations.

Based on localized molecular orbitals, the proposed method reduces large configuration interaction (CI) spaces while maintaining agreement with reference values. Our strategy concentrates the numerical effort on physically pertinent CI-contributions and is to be considered as a tool to tackle large systems including numerous open-shells. To show the efficiency of our method we consider two 4-electron parent systems. First, we illustrate our approach by describing the van der Waals interactions in the (H2)2 system. By systematically including local correlation, dispersion and charge transfer mechanisms, we show that 90% of the reference full CI dissociation energy of the H2 dimer is reproduced using only 3% of the full CI space. Second, the conformational cis/trans rotation barrier of the butadiene molecule is remarkably reproduced (97% of the reference value) with less than 1% of the reference space. This work paves the way to numerical strategies which afford the electronic structure determination of large open-shell systems avoiding the exponential limitation. At the same time, a physical analysis of the contents of the wave function is offered.

The aggressive phenotype: center-level patterns in the utilization of suboptimal kidneys.

Despite the fact that suboptimal kidneys have worse outcomes, differences in waiting times and wait-list mortality have led to variations in the use of these kidneys. It is unknown whether aggressive center-level use of one type of suboptimal graft clusters with aggressive use of other types of suboptimal grafts, and what center characteristics are associated with an overall aggressive phenotype. United Network for Organ Sharing (UNOS) data from 2005 to 2009 for adult kidney transplant recipients was aggregated to the center level. An aggressiveness score was assigned to each center based on usage of suboptimal grafts. Deceased-donor transplant volume correlated with aggressiveness in lower volume, but not higher volume centers. Aggressive centers were mostly found in regions 2 and 9. Aggressiveness was associated with wait-list size (RR 1.69, 95% CI 1.20-2.34, p = 0.002), organ shortage (RR 2.30, 95% CI 1.57-3.37, p < 0.001) and waiting times (RR 1.75, 95% CI 1.20-2.57, p = 0.004). No centers in single-center OPOs were classified as aggressive. In cluster analysis, the most aggressive centers were aggressive in all metrics and vice versa; however, centers with intermediate aggressiveness had phenotypic patterns in their usage of suboptimal kidneys. In conclusion, wait-list size, waiting times, geographic region and OPO competition seem to be driving factors in center-level aggressiveness.

Double excitations in finite systems.

Time-dependent density-functional theory (TDDFT) is widely used in the study of linear response properties of finite systems. However, there are difficulties in properly describing excited states, which have double- and higher-excitation characters, which are particularly important in molecules with an open-shell ground state. These states would be described if the exact TDDFT kernel were used; however, within the adiabatic approximation to the exchange-correlation (xc) kernel, the calculated excitation energies have a strict single-excitation character and are fewer than the real ones. A frequency-dependent xc kernel could create extra poles in the response function, which would describe states with a multiple-excitation character. We introduce a frequency-dependent xc kernel, which can reproduce, within TDDFT, double excitations in finite systems. In order to achieve this, we use the Bethe-Salpeter equation with a dynamically screened Coulomb interaction W(omega), which can describe these excitations, and from this we obtain the xc kernel. Using a two-electron model system, we show that the frequency dependence of W does indeed introduce the double excitations that are instead absent in any static approximation of the electron-hole screening.

A physical model for the longitudinal polarizabilities of polymer chains.

The aim of this work is to provide a physical model to relate the polarizability per unit cell of oligomers to that of their corresponding infinite polymer chains. For this we propose an extrapolation method for the polarizability per unit cell of oligomers by fitting them to a physical model describing the dielectric properties of polymer chains. This physical model is based on the concept of a dielectric needle in which we assume a polymer chain to be well described by a cylindrically shaped nonconducting rod with a radius much smaller than its length. With this model we study in which way the polarizability per unit cell approaches the limit of the infinite chain. We show that within this model the macroscopic contribution of the induced electric field to the macroscopic electric field vanishes in the limit of an infinite polymer chain, i.e., there is no macroscopic screening. The macroscopic electric field becomes equal to the external electric field in this limit. We show that this identification leads to a relation between the polarizability per unit cell and the electric susceptibility of the infinite polymer chain. We test our dielectric needle model on the polarizability per unit cell of oligomers of the hydrogen chain and polyacetylene obtained earlier using time-dependent current-density-functional theory in the adiabatic local-density approximation and with the Vignale-Kohn functional. We also perform calculations using the same theory on truly infinite polymer chains by employing periodic boundary conditions. We show that by extrapolating the oligomer results according to our dielectric needle model we get good agreement with our results from calculations on the corresponding infinite polymer chains.

Ultranonlocality in time-dependent current-density-functional theory: application to conjugated polymers.

We solve the long-standing problem of the large overestimation of the static polarizability of conjugated polymers obtained using the local density approximation within density-functional theory. The local approximation is unable to describe the highly nonlocal exchange and correlation (xc) effects found in these quasi-one-dimensional systems. Time-dependent current-density-functional theory enables a local current description of ultranonlocal xc effects using the Vignale-Kohn functional [G. Vignale and W. Kohn, Phys. Rev. Lett. 77, 2037 (1996)]. Except for the model hydrogen chain, our results are in excellent agreement with the best available correlated methods.

Living The Discipline on a stem cell transplant unit: spiritual care outcomes among bone marrow transplant survivors.

An oncology chaplain details, using The Discipline for Pastoral Care Giving, common themes in the spiritual journeys of stem cell/bone marrow transplantation survivors, and helpful chaplain interventions.